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ABSTRACT 

Astro-H wiW be the first X-ray observatory to employ a high-resolution microcalorimeter, capable of measuring the 
shift and width of individual spectral lines to the precision necessary for estimating the velocity of the diffuse plasma 
in galaxy clusters. This new capability is expected to bring significant progress in understanding the dynamics, and 
therefore the physics, of the intracluster medium. However, because this plasma is optically thin, projection effects will 
be an important complicating factor in interpreting future Tstro-iJ measurements. To study these effects in detail, we 
performed an analysis of the velocity field from simulations of a galaxy cluster experiencing gas sloshing, and generated 
synthetic X-ray spectra, convolved with model Astro-H Soft X-ray Spectrometer (SXS) responses. We find that the 
sloshing motions produce velocity signatures that will be observable by Astro-H in nearby clusters: the shifting of the 
line centroid produced by the fast-moving cold gas underneath the front surface, and line broadening produced by 
the smooth variation of this motion along the line of sight. The line shapes arising from inviscid or strongly viscous 
simulations are very similar, indicating that placing constraints on the gas viscosity from these measurements will 
be difficult. Our spectroscopic analysis demonstrates that, for adequate exposures, Astro-H will be able to recover 
the first two moments of the velocity distribution of these motions accurately, and in some cases multiple velocity 
components may be discerned. The simulations also confirm the importance of accurate treatment of PSF scattering 
in the interpretation of Astro-H/SX.S spectra of cluster plasmas. 

Subject headings: galaxies: clusters: intracluster medium — techniques: spectroscopic — X-rays: 
galaxies: clusters — methods: numerical 


1. INTRODUCTION 


X-ray observatories have yielded a wealth of infor¬ 
mation about the thermodynamic and chemical proper¬ 
ties of the intracluster medium (ICM) of galaxy clus¬ 
ters. However, measuring the kinematics of the ICM has 
been up to now beyond the capability of present instru¬ 
ments. Gas motions in the ICM can be detected by the 
Doppler shifting and broadening of spectral lines: the 
former caused by large-scale bulk motions, and the lat¬ 
ter by turbulence or a complex projection of components 
with different bulk motions along a line of sight. The Re¬ 
flection Grating Spectrometer (RGS) grating on XMM- 


Newton can provide up] 
of spectral lines in cool- 

per limits on Doppler broadening 

core clusters (Sanders et al.|2011: 

ISanders & Fabian|2013| 

1 Pinto et al.|2U15| and references 

therein), but so far no d 

iirect velocity measurements have 


been made, mainly due to the inadequate spectral reso¬ 
lution of e xisting instruments. 

Astro-H (jTakahashi et al.||2014|), a joint JAXA/NASA 
endeavor, will be launched m early 2016, and will be 
the first X-ray observatory capable of detecting mo¬ 
tions in the intracluster medium of galaxy clusters using 
measurements of the shifting and broadening of spec¬ 
tral lines. Astro-H will possess a Soft X-ray Spec¬ 
trometer (SXS) micro-calorimeter with an energy reso¬ 
lution of AE < 7 eV (FWHM) within the energy range 
E ^ 0.3 — 12.0 keV, covering a 3’x3’ field. At the energy 
of the Fe-K(a line, E ^ 6.7 keV, this enables the mea- 
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surement of velocities at resolutions of tens of km s“^. 

Determining the properties of gas motions in clusters 
is important to studies of cluster astrophysics and cos¬ 
mology. Turbulence and bulk motions can serve as a 
transport mechanism for heat , metals, and cosmic rays 
throughout the cluster (e.g., Fujita et al. 2QQ4| Den¬ 


nis fc Ghandran] I2QQ51 Rebusco et al. " [!^udD| [Vazza et 
al. 2Q1Q| Enfilin et ^ 2Qlip ~ Furthermore, these gas 


motions provide non-negligible pressure support against 
gravity, biasing mass estimates based on the assump¬ 
tion of hydrostatic equilibrium, and likely explain dis- 
crepanc ies between hydrostatic and weak lensing-derived 
masses dZhang et a l.l 2Q1Q| [Mahdavi et al.J |2Q13 
Linden et al. I 20141 1 Applegate et al.||2U14|. biu 


von der 
Simulations 


predict that up to ^2U-3U7o of pressure support in even 
some relaxed clusters is due to non-thermal sources, most 


of which will be comprised of gas motions (Evrard et al. 
19961 IRasia et al~1|2QQ6| pN^agai et al.||2QQ7 IFiffaretti ^ 


iqubt IKasia et aI.||zUUbl llNagai et aI.||2UU/t ILmaret 
Valdarnini 120081 Nelson et al. 2014|). Finally, ICM 
bulence is likely a key ingredient tor the origin of 


tur- 

origm ot non- 


therm al phenomena such as radio halos and radio mini¬ 


halos (Brunetti fc L azarian 2007 Donnert et al. 2013 
ZuHone et al. ' 2013p . X-ray instruments with micro- 
calorimeters such as Ast ro-H, Athen^ and the mission 
concept X-ray Surveyor (Weisskopf et al. 2015), are es¬ 
sential to reveal the details of the kinematics of the clus¬ 
ter gas, which will shed light on these questions. 

One class of clusters, those with bright central “cool” 
cores, appear relatively relaxed. However, in these clus¬ 
ters there is one type of gas motion in cool-core galaxy 
clusters that can already be inferred from imaging stud¬ 
ies: “sloshing” motions, evidenced by the presence of 


^ http://www.the-athena-x-ray-observatory.eu/ 
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sharp discontinuities in surface brightness and temper¬ 
ature. In these features, the denser (brighter) side of 
the discontinuity is colder than the lighter (faint er) side, 
hence they have been d ubbed “cold fronts” (see [Marked 


vitch fc Vikhhnin||2QQ7 , for a review). Typically, one or 


several cold fronts may appear in the core region, often 
laid out in a spiral pattern. Simulations have demon¬ 
strated that these features can be produced by gravita- 


ters with subclusters or strong shocks (jChurazov et al. 

2003 

[Tittley & Henriksen|2005[|Ascasibar X. Markevitch 

2006 

ZuHone et al. 2010| 

Roediger et al. 201 ip. The 


are subsonic, with typical velocities of several hundred 
km s~^ {M ^ 0.3 — 0.5) and associated length scales of 
^lOOs of kpc, even up to Mpc in at least a few clusters 
(iSimionescu et ah 2 Q 12 [ Rossetti et al. 2013 Walker et 


al.||2014 ). in the simulations, these gas m otions amplity 


magn etic fields to near-equipytition levels dZuHone et al 


201 i[), gene rate turbulence (jVazza et al.||2012[ |/ 
et al. 2 QT^, and advect and mix metals and ( 


IZuHone 


.. _ _s and entropy 

( [ZuHone et al.||2010l Roediger et al.|[2M2 ). By some es- 


timates, sloshing cold trouts appear m as many as half 


to tw o-thirds of relaxed galaxy clusters (Ghizzardi et al. 


2010 ), indicating that the associated bulk motions are 


common in the cores of massive, non-merging clusters. 
Therefore, it is important to determine if these motions 
can be detected by Astro- and, if so, how they affect 
the shape of the observed lines. 

Previous works have used simulated high-resolution X- 
ray spectra from simulations of galax y clusters to mea¬ 
sure charac teristics of the gas velocity. Inogamov fc Sun^ 


yaev (2003) performed an early analysis of the effects on 
spectral lines from Do ppler shifting and b roadening by 


a turbulent medium. Fujita et al. (2005) showed that 
X-ray spectra with resolution of several eV could mea- 
sure the turbu l ent ga s velocity in a cluster cool core. 


Rebusco et al. (2008) looked at the radial dependence 
of the line width for isotropic and radially directed gas 
motions, proposing to use this radial depend ence as a 
test of th e directionality of the gas velocity. |Shang fc 


Oh| ( 2012 ) demonstrated that mixing model analyses of 


Astro-H spectra could distinguish between different com¬ 
ponents of the velocity distribution , provi ded the expo¬ 


sure ti me was adequate. Nagai et al. (2013) and Ota et al. 
(2015) used adaptive mesh refinement cosmological sim¬ 


ulations and mock Astro-H observations to measure the 
velocity dispersion profiles of relaxed galaxy clusters and 
line shifts an d widths for a mer ging system with multiple 

2013) used clusters produced in 


Bifh et al. 


components. 

a smoothed particle hydrodynamics cosmological simu¬ 
lation and mock Athena observations to investigate the 
velocity structure of the ICM and its impact on the L — T 
relation. 

Our aim in this work is to examine the specific case 
of sloshing motions in a massive galaxy cluster, and “ob¬ 
serve” them with Astro-H. For this, we use cluster merger 
simulations from (jZuHone et al.|2010| hereafter ZMJIO), 
which included runs with varying viscosity. We also em¬ 
ploy recently developed tools for producing synthetic X- 
ray images and spectra from simulations, and produce 
mock observations which are convolved with the appro¬ 
priate instr umental responses to make them as realistic 
as possible ( ZuHone et al.||2Q14 ). This work will be rele¬ 


vant for a number of nearby clusters with indications of 
gas sloshing, including those that already have planned 
Astro-H obs erva tions and others that merit investigation 
(see Section 4.2). 

This paper is organized as follows: in Section we de¬ 
scribe the cluster merger simulations and our m^Ehod for 
producing synthetic X-ray observations. In Sectionwe 
first examine the characteristics of the velocity field Rom 
the simulations, and then generate synthetic X-ray spec¬ 
tra, performing standard spectral analyses to determine 
what aspects of the velocity field may be measured from 
its effect on the spectral lines. We then compare these 
results to those from the simulation. Finally, in Section 
[^we discuss implications of our results and limitations 
of our methods, and in Section we summarize our con¬ 
clusions. We assume a ACDM cosmology with h = 0.71, 
Qrn = 0.27, and Qa = 0.73. 

2 . METHODS 

2 . 1 . Hydrodynamic Simulations of Cluster Mergers 

For our cluster candidate, we examine a minor merger 
simulation from ZMJIO, an off-center collision between 
a large, cool-core cluster and a smaller gasless subclus¬ 
ter, where it is assumed that the smaller cluster has been 
stripped of its gas on a previous passage]^ This configu¬ 
ration produces sloshing cold fronts in the large cluster’s 
core. This simulation w as performed using the FLASH 
code (Dubey et al. 2009), a parallel hydrodynamics/W 
body adaptive mesh refinement (AMR) astrophysical 
simulation code. Full details of the physics employed 
and the initial setup of the simulation can be found in 
ZMJIO, but we provide a short summary here. 

The simulation used FLASH’S standard hydrodynam- 
ics module employing the Pi ecewise-Parabolic Method of 
Colella & Woodward (|1984|) for treatment of the cluster 
plasma, under the assumption of an ideal gas equation 
of state with 7 = 5/3, and a mean molecular weight of 
ji = 0.592, appropriate for an ionized H/He gas with a 
hydrogen mass fraction of X = 0.75. We will examine 
two versions of this simulation: one inviscid, and an¬ 
other with a significant physical viscosity. In th e latter, 
we used an isotrop ic Spitzer dynamic viscosity (Spitzer 

T^[Sarazin|[TMt : 


/i = 0.960 


n{k-QT 


( 1 ) 


rj^bl2 

!2.2 X 10“^®!-— g cm“^ s“\ 

In Ai 


where the temperature T is in Kelvin and the ion 
Coulomb logarithm InAi « 40, appropriate for condi¬ 
tions in the ICM. It is unlikely that the ICM is this vis¬ 
cous, due to the anisotropic n ature of the ion v iscosity in 
a high-/3 magnetized plasma ( Braginskii||1965 ), and also 
because microscale plasma instabilities niay set an upper 
limit on the viscosity tha t is much lower th an expected 
for a collisional plasma (Kunz et al. 2014). However, 
this simulation still serves as a useful test case, since it 
allows us to examine the effects of the sloshing motions 

^ A gas-filled subcluster with the same orbital setup would pro- 
duce a shock front an d turbulence without smooth cold fronts, see 
[ZuHone et al.| ( |2010| ). 
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Fig. 1. — Projections of various quantities along the z-axis of the inviscid simulation, ~1.6 Gyr after core passage. Upper-left: X-ray 
surface brightness (0.5-7 keV band), upper-right: spectroscopic temperature, lower-left: line-of-sight velocity, lower-right: line-of-sight 
velocity dispersion. Square dashed regions indicate the locations of simulated Astro-H SXS pointings (FOV of 3x3), and solid elliptical 
regions correspond to the locations of regions of interest from which we will examine velocity distributions and extract spectra. The cyan 
dot-dashed lines indicate the locations of slices shown in Figures and 


on the spectral lines in the limit that turbulence and in¬ 
stabilities are completely suppressed. The collisionless 
dark matter component of the two clusters is modeled as 
a set of gravitating particles, using an A/'-body module 
which uses the particle-mesh method to map accelera¬ 
tions from the AMR grid to the particle positions. The 
gravitational field due to both the ga s and dark m atter 
is computed using a multigrid solver (Ricker||2008). 

Our merging clusters consist of a large, “main” cluster 
with a mass of 1.25 x 10 ^^ Mq, and a small infalling sub¬ 
cluster, with a mass of 2.5 x 10^^ Mq, for a ma ss ratio of 


5. Our ma in cluster closely resembles A2029 (Vikhlinin 
et al.|[2QQ5 ), a hot, relatively relaxed cluster with slosh- 


mg m the cool core. Initially, the two clusters are set 
up in hydrostatic and virial equilibrium, at a distance of 
3 Mpc away from each other, with an impact parameter 


of 0.5 Mpc on a bound orbit. Each simulation contains 
^10 million dark matter particles and has a finest cell 
size of Ax r\j 5kpc. 

2.2. Simulation of X-ray Photons and Synthetic 
Observations 


To generate X-ray photons for our simulated observa- 
tions, we employ the photon_simulator analysis m odule 
(IZuHone et al.|2QT4 ) from the yt software package (Turk 
et al. 201 l[ ). We~m. 


will outline in brief the procedure for 
generat ing these observation s here, but the reader should 
consult |ZuHone et al.|(|2014|), as well as|Bifii et al.|(|2012|), 
which presents the original implementation of this algo- 
rithm, for more details. 

A relevant 3D region of cells are selected, within which 
the X-ray sp ectrum is compute d for each cell using an 
APEC model ( Smith et al.||2QQl ), using the cell’s density 
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Fig. 2. — Projections of various quantities along the a:-axis of the inviscid simulation, ~1.6 Gyr after core passage. Upper-left: X-ray 
surface brightness (0.5-7 keV band), upper-right: spectroscopic temperature, lower-left: line-of-sight velocity, lower-right: line-of-sight 
velocity dispersion. Square dashed regions indicate the locations of simulated Astro-H SXS pointings (FOV of 3x3), and solid elliptical 
regions correspond to the locations of regions of interest from which we will examine velocity distributions and extract spectra. The cyan 
dot-dashed line indicates the location of the slices shown in Figure 


and temperature, with the effect of thermal broaden¬ 
ing on the emission lines included. We assume a spa¬ 
tially uniform metalli city of Z = 0.3 usin g the 
abundance ratios from Anders & Grevesse (1989). The 
model clusters are situated at the redshift ^^cosmic = 0.05, 
which sets the angular diameter distance of the source. 
The corresponding angular scale is ^ 0.965 kpc/arcsec 
57.88 kpc/arcmin). Initially, a large number of pho¬ 
tons is generated within the 3D region, by assuming very 
large, unrealistic values for the exposure time and effec¬ 
tive area of the instrument. We choose values for these 
parameters such that we obtain roughly ^5-10 x more 
photons than will be eventually used in our synthetic 
observations, allowing us to use these initial photons as 
a statistically representative Monte Carlo sample from 
which to draw subsets for individual exposures. 


Next, we project the photons along the desired line 
of sight, reducing our 3D position space to 2D, and we 
Doppler-shift the photon energies using the line-of-sight 
velocity in the cells from which they originated. In this 
step, we draw a subset of photons from the original sam¬ 
ple, corresponding to a more realistic exposure time (the 
effective area will be determined by the instrumental re¬ 
sponses in the next step). The photon energies are cos- 
mologically redshifted. Lastly, we use t he Tnebingen- 
Boulder ISM absorption model (tbabs, [Wilms et al. 


2000) to model Galactic foreground absorption, assum- 


2 X 10^^ cm 


ing a Galactic column density of Nh = 

To produce synthetic Astro-H observations from our 
simulated event lists, we use The SIMX soft- 

^ http://hea-WWW.harvard.edu/simx/ 
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Fig. 3.— Projections of various quantities along the y-axis of the inviscid simulation, ~1.6 Gyr after core passage. Upper-left: X-ray 
surface brightness (0.5-7 keV band), upper-right: spectroscopic temperature, lower-left: line-of-sight velocity, lower-right: line-of-sight 
velocity dispersion. The square dashed region indicates the location of a simulated Astro-H SXS pointing (FOV of 3x3), and the solid 
elliptical region corresponds to the location of a region of interest from which we wil l ex amine the velocity distribution and extract a 
spectrum. The cyan dot-dashed line indicates the location of the slices shown in Figure 


ware package is capable of simulating the instrumen¬ 
tal response of photon-counting detectors on a num¬ 
ber of current and future X-ray missions, including 
Astro-H. SIMX simulates the predicted effective area, 
vignetting, PSF, detector response, and pileup frac- 
tiori3 of the Astro-H satellite. We use it in con¬ 
junction with our simulated events to produce syn¬ 
thetic SXS spectra. For instrumental responses, we 
use a 5 eV resolution redistribution matrix file (RMF, 
ah_sxs_5ev_20130806 .rmf) and an auxiliary response 
file (ARF, sxt-s_140505_ts02um_intallpxl. arf ), dis¬ 
tributed with SIMX. 

The outputs of SIMX are standard OGIP FITS event 
files which may be viewed and analyzed with standard 
tools (e.g., ds9, XSPEC, FTOOLS, CIAO). SIMX is also ca- 

® for SXS only 


pable of simulating the instrumental background of SXS, 
but we have determined that for our case this background 
emission is an order of magnitude or more smaller than 
the continuum emission in the energy range of interest 
(6.0-7.0 keV, surrounding the Fe-K lines), so for simplic¬ 
ity it is not included in our simulations. We have also 
verified that a typical astro physical backgroun d from un¬ 
resolved AGN sources (see Bantz et al.||2009 , for an ex¬ 
ample) is also significantly smaller than the continuum 
emission for our source, so we have also not included it 
in our simulations. 


3. RESULTS 

As a first step, we will examine projections of the veloc¬ 
ity field obtained from the raw simulation data, to iden¬ 
tify regions which may produce observable consequences 
in our spectral analysis. Second, we will examine the 
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Fig. 4.— Projections of various quantities along the z-axis of the viscous simulation, ~1.6 Gyr after core passage. Upper-left: X-ray 
surface brightness (0.5-7 keV band), upper-right: spectroscopic temperature, lower-left: line-of-sight velocity, lower-right: line-of-sight 
velocity dispersion. Square dashed regions indicate the locations of simulated Astro-H SXS pointings (FOV of 3x3), and solid elliptical 
regions correspond to the locations of regions of in ter est f rom which we will examine velocity distributions. The cyan dot-dashed lines 
indicate the locations of slices shown in Figures and 


properties of the velocity field in detail within specific 
regions. Finally, we will construct and analyze synthetic 
spectra, making connections where appropriate to our 
results from the previous steps. 

3.1. Velocity Moment Maps 

Before constructing our synthetic observations, it is in¬ 
structive to examine the velocity statistics of the raw sim¬ 
ulation data along various lines of sight. This will serve 
as a guide to which regions and projections may yield in¬ 
teresting and observable consequences for Astro-H spec¬ 
tral analysis of our simulations and future observations 
of sloshing cold fronts. 

Throughout this work, we will refer to the cold gas 
component of sloshing cold fronts as situated “below” the 
front surface, and the hot component as the gas “above” 
the front surface, where the sense of “up” and “down” is 


defined with respect to the direction of the local gravita¬ 
tional acceleration, pointing towards the cluster poten¬ 
tial minimum. The sloshing motions have a particular 
geometry which immediately suggests along which lines 
of sight the velocity field would be expected to produce 
the most observable consequences. The mutual orbit of 
the main cluster and subcluster is situated within the 
X — y plane, so most of the gas motion is in the x and y 
directions. Simulations also show that there is a general 
radial expansion of the cold fronts away from t he clus¬ 


ter po tential minimum (Roediger et al. 2011). Keshet 


(2012|) demonstrated that such expanding, spiral flows 
must nave a “bulging” cylindrical or “open barrel” shape 
as seen from lines of sight parallel to the sloshing plane. 

Figures Si and illustrate these characteristics by 
showing projections along the 2 ;, x, and y axes of the in- 
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Fig. 5.— Projections of various quantities along the ar-axis of the viscous simulation, ~1.6 Gyr after core passage. Upper-left: X-ray 
surface brightness (0.5-7 keV band), upper-right: spectroscopic temperature, lower-left: line-of-sight velocity, lower-right: line-of-sight 
velocity dispersion. Square dashed regions indicate the locations of simulated Astro-H SXS pointings (FOV of 3x3), and solid elliptical 
regions correspond to the locations of regions of i nte rest from which we will examine velocity distributions. The cyan dot-dashed line 
indicates the location of the slices shown in Figure pT] 


viscid simulation at an epoch t ^ 1.6 Gyr past the core 
passage. Since the 2 ;-axis projection shows the character¬ 
istic spiral shape of the cold fronts most clearly, we show 
it first. The top panels of these figures show the X-ray 
surface brightness (top left) and gas temperature (top 
righ t), the latter using a “spectroscopic-like” weighting 
Ws\ ( Mazzotta et al.|[2QQ4[ ): 


j(U=y T{r)wsi{r)h' 

^ 2 ^— 3/4 


dr 


Ws\(Xp 


( 2 ) 

( 3 ) 


where h is the unit normal vector defining the line of 
sight direction, and x is the 2-D position vector on the 
sky. The cold fronts appear in all projections, though in 
the X and ^-axis projections there is no evidence of a spi¬ 
ral pattern. However, we will see that these projections 


show the strongest evidence of gas motions from Doppler 
shifting and broadening of spectral lines. 

The bottom left panels of Figures and show the 
line of sight velocity pn (the line shut): 

= J Vnir)wFeir)n-dr (4) 

where in this case the weighting function rcpe is propor¬ 
tional to the emission in the He-like Fe line at ^6.7 keV 
in the cluster rest frame: 

rcpe OC 5Fe(F, Z) (5) 

Throughout this work, the sign convention for line-of- 
sight velocities is such that gas with a positive (negative) 
velocity is moving toward (away from) the observer. 

As expected, we find fairly large (but still subsonic) 
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Fig. 6. — Projections of various quantities along the y-axis of the viscous simulation, ~1.6 Gyr after core passage. Upper-left: X-ray 
surface brightness (0.5-7 keV band), upper-right: spectroscopic temperature, lower-left: line-of-sight velocity, lower-right: line-of-sight 
velocity dispersion. The square dashed region indicates the location of a simulated Astro-H SXS pointing (FOV of 3x3), and the solid 
elliptical region corresponds to the location of a region o f int erest from which we will examine the velocity distribution. The cyan dot-dashed 
line indicates the location of the slices shown in Figure 


line-of-sight velocities, // ^ 300 — 500 km/s, for projec¬ 
tions parallel to the x — y plane of the simulation. The 
regions with the highest velocity shifts are located un¬ 
derneath the cold front surfaces, indicating fast motion 
of the cold gas underneath the fronts. In the z-axis pro¬ 
jection, perpendicular to the sloshing plane, the veloc¬ 
ity shifts are smaller, but fiz can vary by as much as 
^200 km s~^ on small length scales (around 10 kpc), 
indicating the presence of turbulence. 

The bottom right panels of Figures and show 
the line of sight velocity dispersion <7^ line width): 

o'n(x) = j i'^(r)u;Fe(r)n-dr - nl{x) (6) 

Moderate velocity dispersions (a ^ 200—300 km/s) are 
observed in the x and ^-axis projections, mostly spatially 


coincident with the regions of large velocity shifts, again 
underneath the cold front surfaces. In the z-axis projec¬ 
tion, similar velocity dispersions are also observed, both 
in the central core region, and underneath the surface of 
the largest cold front. 

Figures show maps of the same quantities at the 
same epodi as Figures p!p] for the viscous simulation. 


The surface brightness and temperature maps reveal cold 
fronts that are sm ooth and free of instabilities, as seen 


in pr evious works (ZuHone et al. 2010 Roediger et al. 


2013). The maps of line-of-sight velocity and velocity dis- 


persion exhibit a lack of small-scale structure, indicating 
that turbulence is indeed strongly suppressed. However, 
the line shifts of the gas components underneath the cold 
front surfaces are similar to those in the inviscid simula¬ 
tion, and the velocity dispersion in these regions is still 
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Fig. 7. — Characteristics of the velocity field along the z-axis of the inviscid simulation, for region “zl”. Upper panels: slices through the 
a: —z-plane at the center of region “zl”, of temperature (left) and the 2 ;-component of the velocity (right). Black lines indicate the center and 
edges of the elliptical cylinder corresponding to the region in Figure^ Lower-left panel: Phase space plot showing the fraction of emission 
as a function of position and velocity within the cylinder. The black nne indicates the emission-weighted average value. Lower-right panel: 
Effect of plasma motion on a “toy” He-like iron line for the emission with the region. 


fairly significant. 

3.2. Velocity Distribution 
3.2.1. Inviscid Simulation 

The moment maps in the previous section indicate a 
number of interesting locations associated with sloshing 
cold fronts that produce line shifts and broadening that 
will be observable by Astro-H. We will use these maps 
to choose locations for Astro-H pointings, and elliptical 
regions within these pointings from which spectra will 
be extracted. We will take advantage of the velocity 
structure revealed by Figures [T][^ from the simulation to 
decide where to locate the pointings. Of course, we will 
not know a priori in real observations where to point in 
this fashion, but the moment maps show that the regions 
with the most significant line shifting and broadening 


are located just underneath the cold front surfaces, the 
locations of which will be available from higher-resolution 
imaging of the same clusters from Chandra and XMM- 
Newton. 

The chosen pointings and regions for the inviscid sim¬ 
ulation are shown in Figures p!][3l Each elliptical region 
on the maps in 2D defines a c^indrical region in 3D, ex¬ 
tended along the line of sight across the cluster. To get 
a sense of the underlying velocity distribution sampled 
by these regions, we examine the velocity field within 
these cylinders. Figures through show the results 
of this exercise. In the remainder ofThis work, we label 
each region by the axis of projection and the number of 
the pointing, as identified in Figures [T][^ For example, 
pointing 1 on the z-axis projection is labeled “zl”. 

We take Figure as an example, and describe the dif- 
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Fig. 8. — Characteristics of the velocity field along the z-axis of the inviscid simulation, for region “z2”. Upper panels: slices through the 
a: —z-plane at the center of region “z2”, of temperature (left) and the 2 ;-component of the velocity (right). Black lines indicate the center and 
edges of the elliptical cylinder corresponding to the region in Figure^ Lower-left panel: Phase space plot showing the fraction of emission 
as a function of position and velocity within the cylinder. The black nne indicates the emission-weighted average value. Lower-right panel: 
Effect of plasma motion on a “toy” He-like iron line for the emission with the region. 


ferent panels of this figure in detail; the descrption will 
apply to all of the related figures. Figures and [^exam¬ 
ine the 2 ;-couiponent of the gas velocity. The top panels 
show slices through the x — z plane of the temperature 
(left) and 2 ;-component of the velocity (right). The slice 
is taken at the y coordinate of the center of the cylindri¬ 
cal region in Figure (where lines showing the location 
of the slice planes for each projection are shown). The 
center and edges of the cylinder within the slice plane 
are indicated by black dot-dashed and dashed lines, re¬ 
spectively, in Figure The lower-left panel of Figure 
shows a plot of the velocity distribution phase space 
in the cylinder, where the color indicates the amount of 
helium-like iron line emission as a function of both the 
position along the line of sight and the line of sight veloc¬ 
ity. The solid black line indicates the emission-weighted 


average velocity along the length of the cylinder. Finally, 
the bottom-right panel of Figure shows the shape of a 
“toy” helium-like iron line, in the absence of nearby lines 
or continuum, computed from the total emission within 
the cylindrical region. The blue lines show the line shape 
due to thermal broadening without velocity broadening, 
whereas the green lines show the combined effects of ther¬ 
mal broadening and plasma motions. 

In Figure region “zl” is situated within the cold 
fronts enveloping the cluster center. The 2 ;-velocity slice 
plot and phase plot indicate the gas motion is mostly ran¬ 
dom on scales smaller than the core size within this re¬ 
gion, with no significant velocity shift, but a range of gas 
motion between v ^ —400 km/s to ^ +400 km/s. The 
line emission plot confirms this, showing a symmetric 
Gaussian-looking line that has been slightly broadened 
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Fig. 9. — Characteristics of the velocity field along the x-axis of the inviscid simulation. Upper panels: slices through the x — y-plane at 
2 ; = 0, of temperature (left) and the ^-component of the velocity (right). Black lines indicate the center and edges of the elliptical cylinders 
corresponding to the regions in Figure!^ Middle panels: Phase space plots showing the fraction of emission as a function of position and 
velocity within the cylinder. The blacK line indicates the emission-weighted average value. Lower panels: Effect of plasma motion on a 
“toy” He-like iron line for the emission within the regions. 
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Fig. 10.— Characteristics of the velocity field along the y-axis of the inviscid simulation. Upper panels: slices through the x — y-plane at 
2 ; = 0, of temperature (left) and the y-component of the velocity (right). Cyan lines indicate the center and edges of the elliptical cylinder 
corresponding to the region in Figure Lower-left panel: Phase space plot showing the fraction of emission as a function of position and 
velocity within the cylinder. The blacK line indicates the emission-weighted average value. Lower-right panel: Effect of plasma motion on 
a “toy” He-like iron line for the emission with the region. 


by the turbulent velocity field. Figure shows an en¬ 
tirely different velocity field in region “z2”. In this case, 
the elliptical region is centered on the outermost part 
of the spiral to the south, and along the ^-direction the 
symmetry of the sloshing motions results in two oppo¬ 
sitely directed gas flows of approximately Vz ^ 200 km/s 
on either side of the x — ^-plane, showed prominently in 
both the 2 ;-velocity slice plot and the phase plot. Though 
the velocity distribution here is very different from that 
of region “zl”, the shape of the line emission is neverthe¬ 
less nearly identical-both are approximately Gaussian in 
shape. 

In Figure now examining the x-component of the ve¬ 
locity, the slice plots in the top panels are taken through 
the X — y plane of the two subclusters’ mutual orbit, 
and cuts through the center of both the “xl” and “x2” 


cylindrical regions. These regions align with the cold gas 
component of the sloshing cold fronts, which are regions 
of significant bulk motion along the x-axis. The phase 
space plots (middle panels) demonstrate that in both re¬ 
gions there is a substantial portion of gas moving along 
the line of sight in each region, causing the line shifts 
seen in Figure There is also a considerable amount of 
variance in the velocity field, both along the line of sight 
(indicated by the large variation in the black line) and 
in the plane perpendicular to the line of sight along the 
length of the cylinder. For example, at the x = 0 posi¬ 
tion in both regions, there is a An ^ 200 km/s spread in 
velocity from the average value. In the bottom panels of 
Figure the large bulk motions along the line of sight in 
each region create line centroid shifts of An ^ 300 km/s 
(AF^ ~ 7 eV) in either direction, with a fairly significant 
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Fig. 11.— Characteristics of the velocity field along the x-axis of the viscous simulation. Upper panels: slices through the x — y-plane at 
2 ; = 0, of temperature (left) and the ^-component of the velocity (right). Black lines indicate the center and edges of the elliptical cylinders 
corresponding to the regions in Figure]^ Middle panels: Phase space plots showing the fraction of emission as a function of position and 
velocity within the cylinder. The blacK line indicates the emission-weighted average value. Lower panels: Effect of plasma motion on a 
“toy” He-like iron line for the emission within the regions. 
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Fig. 12.— Direct comparison of the velocity-broadened lines (see Figures [7| 11 1 and |18l20t between the inviscid and viscous simulations. 
The energy-dependent emission in each line has been renormalized by its total emission. 


line broadening of FWHM ^ 400 — 500 km/s in both 
cases. Each line also appears slightly skewed in the di¬ 
rection of the velocity shift. These features of the lines 
may be identified with features of the phase plots, which 
both show a broad distribution of gas in velocity space 
centered around a non-zero mean velocity in either di¬ 
rection. 

Lastly, Figure shows the same plots for the single 
region “yl”, for the ^-component of the velocity. The 
slice is taken through the x — y plane at 2 : = 0. Within 
this region, there is a large bulk motion in the direc¬ 
tion. The phase plot shows that within the cold front 
region the bulk of this gas moves with an average y- 
velocity of 1 ; ^ +400 km/s, with a significant spread 
of ^ +200 h800 km/s, seen in the central por¬ 
tion of the panel. A smaller fraction of the X-ray emit¬ 
ting gas is moving in the —y direction with a speed of 
u — 150 km/s, with much smaller variation, seen at a 
distance of ^200 kpc on either side of the center. The 
shape of the emission line in the lower-right panel of Fig¬ 
ure!^ shows evidence of both of these components. 

3.2.2. Viscous Simulation 

We also perform the same examination of the viscous 
simulation. Our results are very similar to those in the 
previous section, so for brevity we only show the results 
for the x-component of the velocity here in Figure pT| and 
refer the reader to Appendix for the figures showing 
the other projections. 

Figure [ 1 TI shows the slice, phase space, and line shape 
plots for the viscous simulation for the x-component of 
the velocity. Overall, due to the damping of turbulence 
and instabilities by viscosity, the velocity distribution 
within each region is smoother, though the large-scale 
motions remain essentially the same. The slice plots are 
almost devoid of small-scale features, and though the 
average value of the velocity along the length of each 
cylindrical region is roughly the same as in Figure]^ (the 
black lines in the phase plots), the emission is not quite 
as spread out in phase space as in the inviscid simulation. 
This manifests itself in a slight narrowing of the velocity 
distribution at any given position along the axis in the 


phase space plots. Nevertheless, the line profiles from 
the viscous simulation are very similar to those from the 
inviscid case. Figure [T^ demonstrates this by comparing 
the line shapes from both simulations directly along all 
three principal axes. This indicates that in both cases 
the gas motion on large scales is not only responsible for 
the shift in the spectral line, but is also largely respon¬ 
sible for the line broadening as well, due to its smooth 
variation along the line of sight. Due to this similarity, 
in the next section we will only generate and fit synthetic 
spectra from the inviscid simulation, for brevity. 

3.3. Synthetic Observations 

We generate our synth etic observations according to 
the procedure in Section |2.2| For each projection, we 
generate 200 synthetic event lists, assuming an exposure 
time of 200 ks. We then perform mock observations at 
the pointing locations shown in Figures extracting 
spectra from the elliptical regions within those pointings 
(the same as those used in the previous section). The 
200 realizations of each region will be used to determine 
the model parameters and their confidence limits. 


3.3.1. Fitting Spectra 

A thermally broadened spectral line with line flux pro¬ 
file Fo{E) in the rest frame of the gas that is Doppler 
broadened to the new line profile F{E) can be repre¬ 
sented by: 


F{E) = 




( 7 ) 


which is essentially a convolution with the velocity dis¬ 
tribution function f{v). We choose to model the ther¬ 
mal spectrum and the velocity distribution and fit these 
models to our synthetic spectra in XSPEC, using the 
bapecj^ thermal plasma model, which accounts for ther¬ 
mal broadening of spectral lines and allows for velocity 
broadening to be accounted for with a single parameter 
a that represents the standard deviation of a Gaussian 


' http://heasarc.gsfc.nasa.gov/xanadu/xspec/manual/ 

XSmodelApec.html 
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Fig. 13.— Spectrum, fitted models, and residuals around the Fe-K lines of the “xl” region. Red curves indicate the total spectral model, 
and in the lower-left panel, the blue and green curves represent the individual model components for Model 3. In the lower-right panel, 
the photon energies were generated assuming thermal broadening only. 


velocity distribution. We constrain the velocity shift ji 
of the spectral lines using the difference of the fitted red- 
shift parameter z of the bapec model to the cosmological 
redshift of the cluster model, /i = —{z — ^^cosmic)^, where 
the sign reflects our convention that positive velocities 
are toward the observer, and 2:cosmic = 0.05. 

We test three different models for the velocity distri¬ 
bution f{v). The first model for /(n), “Model 1”, is a 
single-valued velocity distribution at the shift parameter 

/i: 

f{v) = 5{v — ii) (Model 1) (8) 

This corresponds to a bapec model with the velocity 
broadening parameter set to zero and frozen. This model 
serves as a test of whether the data is consistent with the 
assumption of no velocity broadening, even with Astro- 
iTs improved spectral resolution. The second model, 
“Model 2”, incorporates velocity broadening by thawing 


the a parameter and fitting it along with the others: 

f{v) = G{v; /i, cr^) (Model 2) (9) 

where G is a normalized Gaussian distribution. The third 
model, “Model 3”, is a weighted sum of two “Model 1” 
components: 

f{v) = wi6{v — /ii) + W 2 S{v — 112 ) (Model 3) (10) 

The total Model 3 therefore represents a single¬ 
temperature and single-metallicity plasma which is a 
mixture of two single-valued velocity components (the 
redshift parameters) with different normalizations. In 
XSPEC, this is achieved by adding two bapec models, ty¬ 
ing together the temperature and abundance parameters 
of each component, and allowing the redshift and normal¬ 
ization parameters of each component to vary separately. 
In this model, the velocity broadening parameter of each 
component is frozen at zero. Each of our three models 
also incorporates a tbabs absorption model, where the 



































































































16 


Nh parameter is held fixed at the value of 2 x 10^^ cm“^. 
All other parameters are free to vary, unless noted above. 

We minimize the C-statistic ( Cash|19^ ) to determine 
model parameter values. Though the G’-statistic does not 
have utility as a goodness-of-fit test as the x^-statistic 
does, is works properly with Poisson statistics, and does 
not require rebinning the spectrum to approximate Gaus¬ 
sian statistics. This enables us to use the full spectral 
resolution of Astro-H in our model fits. The spectra are 
fit within a broad energy band of 0.3-10.0 keV. 

We also must ensure that our observations are long 
enough so that our estimates of the line shift and 
width have reasonable statistical accuracy. Ota et ^ 
(2015) showed that ^ 200 counts in the He-like iron 
hue complex are required to achieve 20% accuracy on 
the measurement of the turbulent velocity for values of 
cr ^ 200 km s“^. Table la shows the number of counts in 
this line complex for each of our spectra, all of which eas¬ 
ily exceed this requirement. In real clusters, the metal- 
licity in the core region will be nearly solar, instead of 
the 0.3 assumed in this work, so this requirement will 
be fulfilled for even shorter exposure times in the centers 
of clusters. 

It should also be noted that our procedure does not 
take into account other systematic uncertainties that will 
be important. These include the systematic uncertainty 
on the line shift due to uncertainty in the SXS gain sta¬ 
bility, and the uncertainty on the line width due to uncer¬ 
tainty in the line spread function. The gain uncertainty 
is of particular conce rn. Its required accuracy is 2 eV, 


with a g oal of 1 eV (Mitsuda et al. 2014 Takahashi et 
al.|2014 ). These values would correspond to an accuracy 


on the line shift of Augys ^90(45) km s“^ at 6.7 keV. In 
what follows, we will see that this uncertainty will dom¬ 
inate the statistical uncertainty on the line shift. Since 
we are focused in this work on identifying the systematic 
effects arising from projection and non -Gaussian intrin¬ 
sic lin e shapes, we refer the reader to [Kitayama et al. 


(2014) for a thorough discussion of the instrumental and 
calibration uncertainties. 


effect, for each region we performed a second set of oth¬ 
erwise identical spectral simulations, except that we have 
turned off the effect of vignetting and artificially reduced 
the FWHM of the SXS PSF to 0.01”, to compare to our 
uncorrected spectra. For spectral analysis of real obser¬ 
vations, such an artifical reduction of the PSF scattering 
is obviously unavailable, so the core region and the cold 
front regions will have to be modeled concurre ntly. W e 
show an example of such an analysis in Section [3.3.3 
The value of the reduced C-statistic for each model 
fit is tabulated in Table [^ for spectral simulations with 
and without PSF scatteringFigure shows example 
spectra and fitted models for the “xl” region, within the 
E 6.2 — 6.8 keV band, which surrounds the Fe-K lines 
(at our redshift of ^^cosmic = 0.05). We will describe each 
model in turn, referring to the figure as needed. 


TABLE 2 

C-statistic/d.o.f. for Different Model Fits 


Spectrum 

Model 1 

Model 2 

Model 3 

With PSF Scattering 

xl 

10513/9695 

10291/9694 

10312/9693 

x2 

10051/9695 

9982/9694 

9978/9693 

yi 

10984/9695 

10359/9694 

10430/9693 

zl 

10528/9695 

10162/9694 

10174/9693 

z2 

9963/9695 

9941/9694 

9941/9693 

xl* 

10301/9695 

N/A 

N/A 

Without PSF Scattering 

xl 

10829/9695 

10570/9694 

10587/9693 

x2 

10331/9695 

10258/9694 

10255/9693 

yi 

11137/9695 

10567/9694 

10623/9693 

zl 

10542/9695 

10189/9694 

10205/9693 

z2 

9985/9695 

9958/9694 

9959/9693 


*In this case, velocity broadening and shifting has been turned off. 


TABLE 1 

Total Counts in He-like Fe Line Complex 


Spectrum 

Line Counts 

xl 

874 

x2 

498 

yi 

1355 

zl 

4315 

z2 

421 


Another consideration that must be taken into account 
is the broadness of the SXS PSF, with a width of roughly 
^F. Most of our cold front pointings are offset from the 
cluster center, which has a much higher flux, and inspec¬ 
tion of Figures though shows that the velocity struc¬ 
ture of the core is generally very different from that of the 
cold front regions. A number of photons will be scattered 
by the PSF from the central region into our offset point¬ 
ings, biasing the line shift and width. To account for this 


From Table we see that the fit to Model 1 is gener¬ 
ally poor (except in the “z2” case). The upper-left panel 
of Figure shows the spectrum for region “xl” with its 
fitted Mood 1, showing significant fit residuals around 
the Fe emission lines. For most of the fits to Model 1 
for the various regions, we are able to recover the correct 
temperature, since it is strongly constrained by the con¬ 
tinuum emission, but the metallicity is always underesti¬ 
mated by approximately 10-15%. Taking the “xl” region 
as an example, we find that we can recover the expected 
temperature from the simulation of kTsim = 6.18 keV, 
with kTfit = 6.23lo o6 keV, but the metallicity is under¬ 
estimated, with Zsim = 0.3 Zq and = 0.266io oo 9 ^o- 
We find similar results if we ignore PSF scattering. The 
fitted value of the velocity shift is strongly biased if PSF 
scattering is included, with //fit = -145^12 km s“^, far 

® The XSPEC implementation of the C-statistic is normalized 
such that it provides a y^-like goodness-of-fit value in the limit of 
Gaussian statistics, hence we provide both the value of the statistic 
and the number of degrees of freedom. See http://heasarc.gsfc. 
nasa.gov/xanadu/xspec/manual/XSappendixStatistics.html for 
a more detailed discussion. 
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TABLE 3 

Model 2 Parameters and Simulation Values 


Spectrum 

Tfit 

Fsim 

Zfit 

^sim 

Mfit 

Msim 

c^fit 

^^sim 


(keV) 

(keV) 

(Z©) 

(Z©) 

(km/s) 

(km/s) 

(km/s) 

(km/s) 

With PSF Scattering 

xl 

6 20 +0-05 
D.ZU _Q 05 

6.18 

0 203 +0-009 
u.zyo _o.o09 

0.3 

-156 til 

-221 

203 +11 

197 

x2 

5 41 +0-07 
-0.06 

5.21 

0 202 +0-012 
u.zyz _o.oi2 

0.3 

156 til 

206 

157 tit 

149 

yl 

0 19 +0.04 
-0.04 

6.16 

0 200 +0-007 
u.zyo _o.oo7 

0.3 

231 til 

290 

291 tl§ 

254 

zl 

0 3^^ +0.03 

U.OD _o.02 

6.40 

0 205 +0.004 
U.ZDO _0.004 

0.3 

15 ti 

16 

134 tt 

139 

z2 

c' +0.06 

-0.07 

5.73 

0 202 +0-014 
u.zyz o.o08 

0.3 

_Q +13 
^ -11 

-10 

118 til 

110 

Without PSF Scattering 

xl 

0 19 +0-0^ 
-0.05 

6.18 

0 2Q7 +0-009 
U.zu/ _o.o09 

0.3 

-240 tli 

-221 

203 +>2 

197 

x2 

c; 1 Q +0.05 
-0.06 

5.21 

n 2QQ +0-011 
u.zyy -0.013 

0.3 

224 til 

206 

147 til 

149 

yl 

0 10 +0-0^ 
-0.03 

6.16 

n 2QS +0.009 
u.zyo -0.008 

0.3 

322 til 

290 

266 +1 

254 

zl 

0 +0.03 

0.00 _o.03 

6.40 

0.295 

0.3 

14 t| 

16 

133 tt 

139 

z2 

^ 03 +0-09 
0.00 0.07 

5.73 

0.292 

0.3 

-6 Vs 

-10 

122 tli 

110 


TABLE 4 

Model 3 Parameters and Simulation Values 


Spectrum 

Tflt 

'Aim 

^fit 

Aim 

Ml,fit 

M2,flt 

11^1,fit 

11^2,fit 


(keV) 

(keV) 

(Z©) 

(Z©) 

(km/s) 

(km/s) 



xl 

6 13 +0-05 
-0.05 

6.18 

0.292 t°o °°ol 

0.3 

-372 til 

7 +31 
-41 

0 60 +0-04 
-0.05 

n 40 +0-05 

_0.04 

x2 

5 1 Q +0-05 
-0.06 

5.21 

0 298 +0-011 
o.zyo _o.oi3 

0.3 

318 ti 

9 +48 
^ -43 

0 68 +0-06 
O'-OO -0.09 

n 22 +0-09 
-0.06 

yl 

6 12 +0-05 
-0.03 

6.16 

0 290 +0-008 
u.zyu -0.007 

0.3 

492 ti 

37 til 

0 61 +0-03 
-0.03 

0 39 +0-03 
-0.03 

zl 

6 35 +0-03 

U.OD -0.03 

6.40 

0 293 +0-004 
u.ZDO -0.005 

0.3 

153 tli 

-101 tli 

0 47 +0-03 
'0-^' -0.05 

0 53 +0-05 

U.DO -0.03 

z2 

5 63 +0-09 
0.02 -0.07 

5.73 

0 291 +0-013 
u.zyi o.oio 

0.3 

-123 tv 

87 V°o 

n p; +0-10 
"O-O -0.16 

0 50 +0-16 
O'-OU -0.10 


away from the expected value of /igim = -221 km s“^. 
Without PSF scattering, we find yUfit = -238li| 
a much more accurate value. This indicates that mod¬ 
eling of the bright core component will be necessary to 
accurately measure line shifts of cold fronts. 

As a sanity check, if we turn off the effects of Doppler 
shifting and broadening from gas velocity when creating 
the synthetic spectrum. Model 1 produces a much better 
fit to the data than in the case where these effects were 
included (see the last row of Table and the lower-right 
panel of Figure [T^ despite the fact that we are fitting a 
single-temperature model to a spectrum from gas with a 
range of temperatures. In this case, the metallicity is cor¬ 
rectly recovered with = 0.292+Q;QQg Zq. Since Model 
1 provides a good fit to this spectrum, but not the broad¬ 
ened one, this confirms that Astro-H is sensitive enough 
to detect the effects of Doppler shifting and broadening 
on the emission lines from the subsonic sloshing motions. 

Fitting the synthetic spectra with Model 2 results in 
a significant reduction in the fit statistic (except in the 
“z2” case, and in the case of “x2” only a modest im¬ 
provement is achieved), as seen in the top-right panel 
of Figure Table shows the fitted parameters from 
Model 2, for spectra with PSF scattering (top rows) and 
those without (bottom rows). For nearly all the spectra, 
the fitted temperature T^t agrees with the “real” temper¬ 
ature Tsim, with the exception of the “x2” region, where 


photons scattered from the core by the PSF bias the 
temperature upward. With PSF scattering turned off, 
we recover the correct temperature for “x2”. The metal¬ 
licity parameters are much closer to the correct value of 
^sim = 0-3 Zq, though they are all still biased lower. 

With PSF scattering included in the simulated obser¬ 
vation, the velocity shifts of the “xl”, “x2”, and “yl” 
regions, all within the sloshing plane, are biased in the 
direction of zero velocity from their expected values by 
^60-70 km s“^, due to the mostly unshifted photons from 
the bright core scattering into the field of view. The ve¬ 
locity shifts for the spectra without PSF scattering are 
much closer to the true values, but they still exhibit a 
small bias of ^20-30 km s“^ from the true lin e shift , in 
the direction toward zero velocity. In Section 3.3.2[ we 
will see that this is due to the fact that the true veloc¬ 
ity distributions are not strictly Gaussian. The effect 
of PSF scattering on the velocity broadening parameter 
a is not quite as severe. The fitted line shift parame¬ 
ters for the “zl” and “z2” regions exhibit no bias due 
to PSF scattering, but this is expected, since the “zl” 
region is centered on the core itself, and the “z2” region 
has nearly the same line shift, due to the symmetry of 
the sloshing motions. In all cases, however, the velocity 
bias from PSF scattering is comparable to or even less 
than that expected from instrumental uncertainties, if 
^'^sys ^ 45(90) km s~^ as noted above. 
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Fig. 14.— Velocity distributions for the five regions from the inviscid simulation with the predicted velocity distributions from Model 2 
and, in two cases, Model 4. Vertical arrows indicate the position of the line shift associated with the (5-function component in Model 4, 
and the error bars on the arrows give the 1-a error on the line shift. 


Fitting the spectra with Model 3 also provides a good 
fit in most cases (see the bottom-left panel of Figure 
13 for the fit for region “xl”). The reduction in the 
U^statistic is approximately the same as in Model 2. 
Though it might not be expected that such a crude ve¬ 
locity distribution function (Equation [Iq|) would provide 
a good fit to complex velocity distributions, this is made 
possible because of the blending together of the individ¬ 
ual line components of complexes (such as Fe near 7 keV 
in the rest frame) due to the finite spectral resolution of 
SXS (^5 eV) and line broadening. Table 0 shows the 
fitted parameters from Model 3, with no PSF scattering. 
In general, the temperature parameters in each spectrum 
are correctly recovered, and the two velocity components 
/ii and JU 2 have comparable normalization, with a typ¬ 
ical difference of ^250-400 km s“^. Though the model 
provides a good fit, its utility as a physical description of 
the velocity distribution is limited. 


3.3.2. Predicted Velocity Distributions 

How well do our models reproduce the actual velocity 
distributions within the regions from which the spectra 
are extracted? Figure M shows the predicted velocity 
distribution function f{y) from Model 2 (green curves) 
for our five regions, compared to the actual f{v) for the 
same regions (black curves). For the “xl”, “x2”, and 
“yl” regions, the velocity distributions are complex, and 
exhibit deviations from Gaussianity, though the Model 
2 fits do approximately capture the width of the distri¬ 
bution. This explains why the Model 2 line shifts from 


Table (without PSF scattering) do not always agree 
precisSy with the line shift measured from the simula¬ 
tion (though this disagreement is small compared to the 
expected systematic error due to uncertainty in the line 
spread function and gain); the mean of the best-fit Gaus¬ 
sian component and the mean of the true distribution 
are not the same if the underlying distribution is non- 
Gaussian. The velocity distributions for the “zl” and 
“z2” regions are well-fit by Gaussian distributions, and 
as a result both the line shift and width are correctly 
recovered by Model 2. 


Figure m indicates that the velocity distribution of 
cold fr onts may be bette r modeled by multiple compo¬ 
nents. iShang & Oh (2012) used a mixing-model approach 


to fit velocity distribution models to Doppler-broadened 
spectral lines. They found that they were able to fit ve¬ 
locity distributions from a variety of simulated velocity 
fields (including those from core gas sloshing) using a 
sum of Gaussian models. They also found that for the 
Astro-H spatial and spectral resolution, typically only 
tw o Gaussian compone nts can be constrained. Gompared 
to IShang & Oh (2012), we are limited by statistics. For 


each of their spectra, they assumed 10^ counts in the sin¬ 
gle He-like iron line at 6.7 keV, corresponding to roughly 
a megasecond exposure of the entire SXS field of view 
pointed at the cluster core for a number of nearby clus¬ 
ters (see their Table 2). As noted above, our observa¬ 
tions of 200 ks exposure typically have ^10^ counts in 
the entire He-like iron line complex (see Table [^. Also, 
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Fig. 15. — Cumulative velocity distributions for the five regions from the inviscid simulation with the predicted cumulative velocity 
distributions from Model 3. 


we extract spectra from smaller regions underneath cold 
front surfaces, which are fainter than the cluster core. 

Though we do not employ a mixing-model approach, 
we can fit a sum of two Gaussian models by extending our 
Model 3 (Equation [Tq|); we thaw the velocity broadening 
parameters of the separate bapec components, allowing 
them to vary in the fit. Along with the temperature and 
metallicity parameters, we have six additional parame¬ 
ters: a normalization, line centroid, and line width for 
each Gaussian velocity component. In all cases, we are 
unable to constrain all of these parameters uniquely with 
our simulated spectra. However, we do find that we can 
constrain one velocity broadening parameter. This de¬ 
fines our “Model 4”, a simple generalization of Model 
3: 

f(v)=wiS(v-fj,i)+W 2 G(v;/j, 2 ,cr 2 ) (Model 4) (11) 

Predicted velocity distributions from Model 4 fits to the 
“x2” and “yl” spectra are shown in Figure as red 
curves with arrows showing the positions oFthe line 
shift parameter along with its 1-a error bar. We are 
only able to constrain parameters for these particular re¬ 
gions, as these cases have the clearest evidence for two 
well-separated velocity components. In these two cases. 
Model 4 resolves both components well. 

To compare Model 3 to the actual velocity distribution, 
it is more instructive to examine the cumulative distri¬ 
bution function (CDF) of the velocity, F{v). For Model 


3, the CDF of the velocity distribution is 

F{v) = wiH{v — 111 ) + W 2 H{v — 112 ) (Model 3) (12) 

where H is the Heaviside step function. Figure \TE\ shows 
F{v) for the five different regions, with the Model 3 pre¬ 
diction overlaid. From this figure, we see that Model 3 
represents a crude, piecewise-constant approximation to 
the cumulative distribution of the velocity, which may 
provide a rough sense of the location and relative impor¬ 
tance of the two most dominant velocity components, 
but in general a Gaussian velocity model or (if possi¬ 
ble) a sum of two models is to be preferred, as they will 
provide a more accurate physical representation of the 
velocity field for a comparable number of parameters. 


3.3.3. Accounting for PSF Scattering Via Joint Modeling 

As mentioned above, in SXS observations of cold fronts 
a non-negligible number of photons may be scattered into 
the field of view from the nearby core region. There¬ 
fore, to achieve accurate measurements of model param¬ 
eters, it will be necessary to model the spectrum of a 
cold front region simultaneously with the spectrum of the 
cluster core. Examples of such a nalyses can be found in 
the Astro-H dustei white paper (Kitayama et al.||20I4), 
where Chandra images were used in conjunction with 
simulated spectra to estimate the effects of PSF scat¬ 
tering. In this section, we perform such a joint analysis 
using the “xl” and “x2” regions as examples. 
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Fig. 16.— Joint fits to cold front and core region spectra. Left panel: Projections of X-ray surface brightness and line-of-sight velocity 
along the x-axis, with the regions for spectral extraction overlaid. The “core” region has been added. Right panel: Example spectra from 
the “xl” and “core” regions with fitted models to both the “core” and “xl” components. 


TABLE 5 

Accounting for PSF Scattering: Joint Fit Parameters and Simulation Values 


Spectrum 

Tfit 

Uim 

Zfit 

^sim 

/^fit 

/^sim 


©sim 


(keV) 

(keV) 

(Z©) 

(Z©) 

(km/s) 

(km/s) 

(km/s) 

(km/s) 

xl 

99 +0.07 
-0.07 

6.18 

0.293 

0.3 

-192 ±\l 

-221 

196 ili 

197 

x2 

c; Qf) +0.06 
D.OU o.09 

5.21 

0.294 

0.3 

176 ±1° 

206 

161 ii 

149 


The left panel of Figure 16 shows maps of the X- 
ray surface brightness and the line-of-sight velocity, with 
the “xl” and “x2” regions chosen for spectral extraction 
(which are the same as before), with an additional “core” 
region that contains the bright cluster core. We then 
perform separate simulations including only the photons 
originating from this core region that are scattered into 
regions “xl” and “x2”. We create spectra from each re¬ 
gion which contain only photons originating in the core. 
We find that the photon flux from the core scattered into 
the “xl” and “x2” regions is about ^20% of the total flux 
of each region. We performed a joint fit in XSPEC of 
the two spectra using a sum of two bapec models, where 
one component models the photons from the core region 
spectrum and the other models those from the cold front 
spectrum. The right panel of Figure pT] shows example 
spectra for the “xl” region and the fitted models. Table 
shows the fitted parameters for the “xl” region result¬ 
ing from the joint fit, where we only show the parameters 
for the cold front model component. The joint fit results 
in an improvement of the estimation of the line shift of 
the cold front component in both regions. 

This effect may be slightly larger than our simulations 
indicate. In reality, clusters have central metallicities 
of ^ Zq or higher, and the metallicity decreases with 
radius. In our simulations we have assumed a spatially 
uniform Z = 0.3 z©. Depending on the difference in 
metallicity between the core and the cold front for a given 


cluster, the PSF-scattered flux from the core into the cold 
front region may be somewhat larger. 

4. DISCUSSION 

4.1. Implications for Estimating the I CM Viscosity 

Our results show that the line shapes produced by the 
sloshing motions are very similiar, regardless of whether 
or no t the ICM is inviscid or very viscous (see Section 
3.2.2 and Figure 12). This is not unexpected, as the 
largest velocities along the line of sight will be associ¬ 
ated with the driving scale of the sloshing (~100 kpc), 
which is much larger than th e viscous dissipation scale. 
Inogamov & Sunyaev (2003) used a study of ICM tur- 
bulence to show that the last gas motions at the driving 
scale will be mostly responsible for the shape of the spec¬ 
tral lines. Sloshing represents an analogous situation, 
producing smaller-scale turbulent motions with smaller 
velocities which will ha ve less of an effect on the shape 
of the spectral line (see Vazza et al.||2012 ZuHone et al. 

Ewen in nearby clusters, the viscous dissipation scale 
will be unresolved by As tro-H, unless this scal e is implau¬ 
sibly large. For example, ZuHone et al. (2015) showed us¬ 
ing simulations of turbulent velocity flelds in a Coma-like 
cluster that SXS’s spatial resolution will be too coarse to 
constrain the dissipation scale of the turbulent cascade. 
Definitive constraints on viscosity from gas motions in 
clusters will likely have to wait for missions with similar 
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Fig. 17.— Chandra images of nearby clusters with cold fronts. 3’x3’ squares are shown to indicate positions of SXS pointing s which could 

potentiall y yield interesting results based on the findings of this work. All panels except the top-left (A2029) are taken from |Kitayama et| 

- 


spectral resolution and better spatial resolution, such as 
Athena and X-ray Surveyor. 

4.2. Applieation to Possible Astro-H Targets: Where 
Should We Point? 

The results from our simulation analysis of sloshing 
cold fronts may be relevant for a number of nearby clus¬ 
ters with similar features. The cold gas component of 
sloshing, located underneath the front surface, is likely 
to be a region where emission lines will be significantly 
broadened. If our line of sight is within the plane of the 
sloshing motions, we also expect to measure a significant 
line shift. Here we identify a few nearby systems with 
cold fronts which could be observed by Astro-H^ to de¬ 
termine if the line shifting and broadening predicted in 
this work can be detected. 

Two nearby clusters which have clear-cut examples of 
sloshing cold fronts are Abell 2029 and Abell 2319. Abell 
2029 is at a redshift of 0.07728, and in this case the cold 
fronts can be entirely confined within a single SXS point¬ 


ing of 3’x3’ (see the top-left panel of Figure [Tt]). For this 
cluster, the SXS pointing should be placed at the cluster 
center, and the different portions of the cold spiral can 
be resolved at a scale of ^1’, though accurate modeling 
of the effects of the PSF will be required. In Abell 2319, 
at a redshift of 0.0557, the cold fronts are larger in the 
sky plane, and can be split across a few SXS pointings 
(an example pointing is shown in the top-right panel of 
Figure [^. In any case, each pointing will need to have 
long enough exposure to ensure that enough counts will 
be within the Fe-K co mplexes to mea sure the line shift 
and width accurately (Ota et al.||2015|). 

Another interesting system is Abell 3667, which pos¬ 
sesses a large cold front to the southeast (see the bottom- 
left panel of Figure 17). Kinematic measurements in 
the X-ray band have the potential in this case to distin¬ 
guish between different merging scenarios for this cluster. 
Abell 3667 is assumed to be a n ongoing major merger 


nearly in the plane of the sky (Rottgering et al. 1997 
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Vikhlimn et al. pOOllPohnston-Hollitt et al.|2QQ8[|Owers 

et al.||2009| IDatta et al.||2014p. |Vikhliiiin et al. (|200ip 
used hydrodynamic arguments to estimate the velocity 
of the cold front to be transsonic or supersonic, support¬ 
ing this in terpretation, though a recent analysis [Datta et] 
al. (2014) suggests the velocity is subsonic. However, if 
there is a significant line-of-sight velocity shift of the gas 
underneath the cold front surface (the pointing outlined 
in red with a dashed line through the center in Figure 
17), this may indicate that we are observing a sloshing 
cold front with our line of sight at least partially aligned 
with the plane of the gas motions, or at least that the 
merger is not strictly in the plane of the sky. Simula¬ 
tions indicate that sloshing cold fronts can also appear 


in major mergers with large impact parameters (Ricker 
fc Sarazin||2QQl Poole et al.||2QQ6[ [ZuHone||2Qll ). 


4'inally, an extensive Astro-h study is already planned 
for the Perseus Cluster {z = 0.0179), the brightest X-ray 
cluster in the sky. Due to its proximity of ^68 Mpc, 
Astro-H will be able to map the core of Perseus with 
excellent resolution, with F ~ 20 kpc. Using a combina- 
tion of ROSAT, XMM-N ewton, and Suzaku observations, 
Simionescu et al. (2012) demonstrated evidence of large- 
scale sloshing motions extending from the cluster core 
to the virial radius. The bottom-right panel of Figure pT| 
shows the residual image of the cl uster core, with planned 
Astro-H/SXS pointings overlaid (jKitayama et al.||2014p, 
in which the large-scale sloshing spiral can be readily 
seen. It appears that the plane of the sloshing motions 
is oriented nearly perpendicular to our line of sight, so 
on the basis of this work we expect to find no significant 
line shifts associated with the cold gas underneath the 
front surfaces, but may find significant line broadening 
from the expansion of the cold fronts in our line of sight. 
However, Perseus is also a site of significant AGN activ¬ 
ity, which is likely driving significant turbulence, which 
will also produce significant line broadening. Teasing the 
effects of sloshing and AGN-driven turbulence apart will 
be an interesting challenge for future simulations. 

5. SUMMARY 


We have carried out a detailed investigation of gas mo¬ 
tions in simulations of a relaxed, cool-core cluster, specif¬ 
ically with regard to how these motions may be observed 
by the Astro-H X-ray observatory. Our main results are 
as follows: 


• Astro-H wiW be able to detect the large-scale bulk 
motions associated with sloshing gas in cool-core 
clusters. These motions will produce significant 
shifting and broadening of spectral lines that will 
be easily measured by Astro-Hs high-resolution 
calorimeter. The regions with the most detectable 
velocities are located underneath the cold front sur¬ 
faces, associated with the cold, dense phase. 


• The line shifts produced by cold fronts will be most 
observable when our line of sight is directly within 
the plane of the sloshing motions. This is also the 
direction where the evidence of sloshing from imag¬ 
ing studies will be the least obvious; cold fronts will 
still be seen, but without an obvious spiral pattern. 
Line widths seen in this orientation will be pro¬ 
duced predominately by the variation of the slosh¬ 
ing motion along the line of sight. In cases where 


a clear sloshing spiral is seen, our line of sight is 
nearly perpendicular to the sloshing plane, and line 
shifts from the sloshing motions should be insignif¬ 
icant. However, significant line broadening may 
be produced by sloshing-driven turbulence or by 
the radial expansion of the cold fronts perpendic¬ 
ular to the sloshing plane. In either case, sloshing 
motions provide evidence of potentially significant 
non-thermal pressure support. 

• Sloshing motions produce line shapes that are well- 
represented by one or more Gaussian components. 
One-component Gaussian velocity models provide 
good fits to our synthetic spectra, regardless of the 
underlying velocity distribution. However, depend¬ 
ing on the degree of non-Gaussianity of the under¬ 
lying distribution, line shifts and widths under this 
assumption may be biased. For long exposures, 
two-component models may be used to fit complex 
velocity distributions with well-separated compo¬ 
nents. 


• The line shapes produced in our inviscid and highly 
viscous simulations are very similar, and for practi¬ 
cal purposes will be indistinguishable. This reflects 
the fact that the medium-to-large spatial scale mo¬ 
tions have the strongest effect on the line shape, 
which are relatively unaffected by viscous dissipa¬ 
tion, even if the IGM is very viscous. For this rea¬ 
son, constraints on viscosity from line shapes pro¬ 
duced by sloshing motions are likely to be limited. 

• Measured line shifts and widths from sloshing cold 
fronts will be affected by systematic effects, such 
as uncertainty in the line spread function and gain. 
In particular, the systematic error on the line shift 
is likely to dominate over other sources of error dis¬ 
cussed in this work, including statistical uncertain¬ 
ties, errors arising from fitting Gaussian models to 
non-Gaussian velocity distributions, and PSF scat¬ 
tering from other regions. 


• Though a smaller effect, measured line shifts and 
widths will also be affected by PSF scattering of 
photons from the bright cluster core. In our case, 
we find that the photons scattered from the nearby 
core into the field of view centered on the sloshing 
cold fronts can bias line shifts and widths by sev¬ 
eral tens of km/s, though the bias on line widths is 
typically smaller. This indicates that spectra from 
cold front regions in clusters should be modeled in 
combination with that from the core regi on to mit¬ 
igat e this effect, as we sho wed in Section 3.3.3 (see 
also Kitayama et al.|[2Q14 ). 


There are a number of unresolved questions that our 
work does not address, and a number of avenues for fu¬ 
ture investigation are suggested by our results. 


• Though our idealized setup provides the oppor¬ 
tunity to determine what effects sloshing motions 
have on spectral lines, in real clusters the situa¬ 
tion is more complicated. For example, the domi¬ 
nant driver of all gas motion, including turbulence, 
in our simulation is the sloshing motions them¬ 
selves. However, even in relaxed systems, there 
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will be other sources exciting plasma motions, such 
as AGN and substructure, so the gas motion pro¬ 
duced in our simulations is likely a lower limit on 
what can be expected. Analysis of a more com¬ 
plicated system, such as a cluster extracted from a 
cosmological simulation, would provide a fuller pic¬ 
ture of the combined effect of turbulence and bulk 
motions on observations of spectral lines. 

• For this work, we chose to analyze a single epoch 
shortly after the core passage of the subcluster, 
which is late enough so that cold fronts have had 
time to develop, but early enough so that they are 
still bright and prominent. A detailed analysis of 
the later epochs of our simulation are needed to 
determine if our conclusions hold as the cold fronts 
expand. We have had a cursory look at the data 
from later times, which suggests that similar results 
will indeed hold. 


• We predict that line shapes produced by slosh¬ 
ing may have a detectable non-Gaussianity. This 
may be modeled with a superpositi o n of Gaus- 
sians, as suggested by Shang & Oh (2012). Al¬ 
ternatively, new models for spectral fitting could 
be developed that incorporate physically motivated 
non-Gaussian velocity distributions. However, our 
ability to constrain these models will be primarily 
limited by the SXS gain uncertainty. 


• A separate analysis of our data using model re¬ 
sponses for the upcoming Athena and X-ray Sur¬ 
veyor missions would be useful to characterize the 
advantages of larger effective area and smaller PSF. 
This is also possible using the SIMX package. 


• It is not clear to what extent the results of our anal¬ 
ysis would apply to cold fronts produced in major 
mergers. In these systems, shocks are also present, 
the bulk motions are significantly larger, and the 
turbulence produced can be significant. A similar 
analysis of a major merger simulation is necessary. 
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APPENDIX 


A. SPECTRAL MODELING VERIEICATION TEST 

We performed a verification test of our procedure for generating and fitting synthetic Aspectra. In this test, 
we ensure we can recover the input plasma temperature and metallicity, as well as the shift and width of spectral lines, 
for a simplified cluster model. Eor our test, we set up a n isothermal, kT = 6 keV cluster with a density profile given 
by a /3-model (Cavaliere & Eusco-Eemiano 1976 1978), with a core radius Vc = 50 kpc, = 1 , and a core electron 
density ric = 0.035 cm”*^. We then add a velocity field to the cluster model using a single Gaussian random field with 
jiiz = 400 km/s and az = 400 km/s. We perform 200 realizations of this velocity field, and from each one we compute 
a realization of the spectrum taken from an entire SXS pointing located at the center of the model cluster. We use 
the 200 different spectral realizations to compute the mean and l-cr errors for each parameter. 

In the case of the normalization parameter 77 , we must take into account that a fraction of the photons within the 
region covered by the SXS pointing will be scattered out of the region due to PSE scattering and vignetting effects. 
To quantify this effect, we perform the same simulation of the synthetic observation with these effects turned off. 
We determine that ^ 7% of the photons have been removed from the region by these effects, and have adjusted the 
expected value of 77 accordingly. The results of the spectral fitting test are shown in Table We find that we are able 
to recover the values of all of the parameters within the l-cr errors. 


TABLE 6 

Results of Spectral-Fitting Verification Test 


Parameter True Value Fitted Value 


kT (keV) 

6 

5 07 +0-03 
-0.04 

z (Zo) 

0.3 

n 900 + 0.006 
u.zyy -0.005 

liz (km/s) 

400 

420 

az (km/s) 

400 

398 ±11 

r] (10“^ cm“^) 

1.522 

1 con +0.003 
J-.OZU _0.004 


B. VISCOUS SIMULATION EIGURES 

In what follows, we present the rest of the slice, phase space, and line shape plots for the viscous simulation. We 
refer back to Eignres[4][^for the locations of the spectral extracting regions and slice planes. 
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Fig. 18. — Characteristics of the velocity field along the 2 ;-axis of the viscous simulation, for region “zl”. Upper panels: slices through the 
a: —z-plane at the center of region “zl”, of temperature (left) and the 2 ;-component of the velocity (right). Black lines indicate the center and 
edges of the elliptical cylinder corresponding to the region in Figure^ Lower-left panel: Phase space plot showing the fraction of emission 
as a function of position and velocity within the cylinder. The black nne indicates the emission-weighted average value. Lower-right panel: 
Effect of plasma motion on a “toy” He-like iron line for the emission with the region. 
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Fig. 19. — Characteristics of the velocity field along the 2;-axis of the viscous simulation, for region “z2”. Upper panels: slices through the 
a: —z-plane at the center of region “z2”, of temperature (left) and the 2;-component of the velocity (right). Black lines indicate the center and 
edges of the elliptical cylinder corresponding to the region in Figure^ Lower-left panel: Phase space plot showing the fraction of emission 
as a function of position and velocity within the cylinder. The black nne indicates the emission-weighted average value. Lower-right panel: 
Effect of plasma motion on a “toy” He-like iron line for the emission with the region. 
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Fig. 20. — Characteristics of the velocity field along the y-axis of the viscous simulation. Upper panels: slices through the x — y-plane at 
2 ; = 0, of temperature (left) and the y-component of the velocity (right). Cyan lines indicate the center and edges of the elliptical cylinder 
corresponding to the region in Figure Lower-left panel: Phase space plot showing the fraction of emission as a function of position and 
velocity within the cylinder. The blacK line indicates the emission-weighted average value. Lower-right panel: Effect of plasma motion on 
a “toy” He-like iron line for the emission with the region. 












